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Abstract: We develop a theory for pricing non-diversifiable mortality risk in an incomplete 
market. We do this by assuming that the company issuing a mortality-contingent claim 
requires compensation for this risk in the form of a pre-specified instantaneous Sharpe ratio. 
We prove that our ensuing valuation formula satisfies a number of desirable properties. For 
example, we show that it is subadditive in the number of contracts sold. A key result is 
that if the hazard rate is stochastic, then the risk-adjusted survival probability is greater 
than the physical survival probability, even as the number of contracts approaches infinity. 
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1. Introduction 

Insurance markets are incomplete for several reasons, and in this paper we focus on 
the incompleteness arising from two sources: (1) The timing of insurance payments is 
generally determined by a jump process - such as, when individuals die or when cars 
crash, and (2) insurers cannot buy and sell insurance contracts frictionlessly, if at all. 
Actuaries traditionally assume that one can eliminate the uncertainty associated with the 
jump process by selling a large number of insurance contracts. By invoking the law of large 
numbers, actuaries thereby replace the random occurrence of insurance payments with a 
deterministic schedule. 

This assumption is the reason that Brennan and Schwartz (1976), or more recently 
Boyle and Hardy (2003), were able to extend the Black-Scholes formula to price a derivative 
instrument with a random maturity date. They assumed that the issuer sold a sufficiently 
large number of policies so that the issuer only had to contend with a deterministic schedule 
of maturing derivatives. This simplifying technique permeates most of the literature that 
combines finance and insurance. While each individual's payment occurs at a random time, 
one assumes that a large portfolio becomes deterministic and, thus, safe from a mortality 
point of view. 

Using the language of modern portfolio theory, under this assumption, the idiosyn- 
cractic risk - that is, the standard deviation per contract - goes to zero if the insurer sells 
enough contracts, and so the risk is diversifiable. However, because the insurer can only 
sell a finite number of insurance policies, it is impossible for the insurer to eliminate the 
risk that the experience will differ from what is expected. The risk associated with selling 
a finite number of insurance contracts is what we call the finite portfolio risk. 

In addition to the finite portfolio risk, there is a risk arising from the fact that the 
parameters underlying the claim occurrence process are random themselves. For example, 
for an insurance payment that is contingent on the survival or death of an individual, 
one models the claim occurrence process as a Poisson process with intensity A, also called 
the hazard rate. When the individual dies, then the process jumps and is "killed" at 
that time. Now, if A is a deterministic function of time, then the risk associated with 
selling a finite number of insurance contracts, each of whose claim occurrence is dictated 
by independent and identically distributed killed Poisson processes, is the finite portfolio 
risk. If A is stochastic, then there is an additional risk that we call the stochastic mortality 
risk, a special case of stochastic parameter risk. Even as the insurer sells an arbitrarily 
large number of contracts, the systematic stochastic mortality risk remains because of the 
second source of incompleteness, namely, the inability to buy and sell insurance contracts 
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frictionlessly. 

One could also define similar risks associated with the size of insurance payments. 
However, for concreteness in this paper, we focus on the rate of claim occurrence and 
leave considerations concerning the size of the payments to future research. In particular, 
we determine how to value the risks associated with a specific contract called a pure 
endowment. A pure endowment pays $1 to an individual at a given time T if the individual 
survives to that time. Pure endowments are important in that they are the building blocks 
of life annuities. They are also analogous to defaultable zero-coupon bonds because of the 
uncertainty of the payment at time T. 

We emphasize that stochastic parameter risk is different from model specification 
risk. For stochastic parameter risk, we assume that we have specified the model precisely, 
including the stochastic process that the model parameter follows. However, if one cannot 
buy and sell insurance contracts frictionlessly, one cannot eliminate the systematic risk 
inherent in this stochastic process. On the other hand, if an insurer were to issue a pure 
endowment contract to infinitely many individuals, each of whose independent time-of- 
death follows the same stochastic hazard rate, then the insurer could completely hedge 
the stochastic mortality risk by selling term life insurance to infinitely many individuals, 
each of whose independent mortality follows the same stochastic hazard rate as the group 
buying the pure endowment. This remark is the large-scale version of the statement that 
if one could sell both a pure endowment and a term life insurance contract to the same 
person (with the same maturity date and same mortality-contingent payment) , then there 
would be no risk. In this case, the price for the combined contract would be the price for 
a default-free zero-coupon bond. 

In this paper, we consider the case in which the issuer of a pure endowment cannot sell 
term life insurance to hedge the stochastic mortality risk. (In future work, we allow the 
issuer to hedge the risk partially by selling life insurance to an individual whose stochastic 
mortality is correlated with that of the purchaser of the pure endowment.) We argue that 
survival probabilities are uncertain and that this uncertainty is correlated across individ- 
uals in a population - mostly due to medical breakthroughs or environmental factors that 
affect the entire population. For example, if there is a positive probability that medical 
science will find a cure for cancer during the next thirty years, this will influence aggre- 
gate mortality patterns. The uncertainty regarding the evolution of survival probabilities 
induces a mortality dependence that cannot be diversified by selling more contracts. This 
risk induces a mortality risk premium that should be priced by the market and whose 
magnitude depends on a representative investor's risk aversion or demanded compensation 
for risk. 
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We are not the first to recognize that mortality rates themselves should be viewed as 
stochastic. Biffis (2005), Schrager (2005), Dahl (2004), as well as Milevsky and Promislow 
(2001) and Soininen (1995), used and calibrated diffusion processes to model the force 
of mortality. Also, demographers and actuaries, such as Lee and Carter (1992), Olivieri 
(2001), and DiLorenzo and Sibillo (2003), developed methods for projecting mortality. 
While some practitioners, such as Smith, Moran, and Walczak (2003), and academics, 
such as Cox and Lin (2004) and Cairns, Blake, and Dowd (2004), implicitly recognize that 
mortality risk is being priced by the market, they did not derive the actual value of this 
risk from first principles. In related work, Blanchett-Scalliet, El Karoui, and Martellini 
(2005) value assets that mature at a random time by using the principle of no arbitrage; 
the resulting pricing rule is, therefore, linear. However, for insurance markets, one can- 
not assert that no arbitrage holds, so we take a different approach to valuing insurance 
contracts. 

We value a pure endowment by assuming that the insurance company issuing the 
contract is compensated for risk via the so-called instantaneous Sharpe ratio of a suitably- 
defined portfolio. Specifically, we assume that the insurance company picks a target ratio of 
expected excess return to standard deviation, denoted by a, and then determines a price for 
a pure endowment that yields this pre-determined a for the corresponding portfolio. Our 
results do not depend on using a particular diffusion for the hazard rate. In future research, 
we plan to calibrate this model to a generalized mean-reverting process for mortality rates, 
similar to the work by Norberg (2004) in the context of interest rates. 

Our methodology recovers a number of results that one expects within the context 
of insurance, but we also obtain new insights into the breakdown of traditional insurance 
pricing. For example, we prove that if the hazard rate is deterministic, then as the number 
of contracts approaches infinity, the price of a pure endowment collapses to the discounted 
expected payment using the physical probability measure to value the mortality risk, re- 
gardless of the target value of the Sharpe ratio. In other words, if the stochastic mortality 
risk is not present, then the price for a large portfolio reflects this and reduces to the 
"usual" expected value pricing rule in the limit. 

A key result of the paper is that if the hazard rate is stochastic, then the financial 
value of the pure endowment is greater than the above-mentioned discounted expected 
payment, even as the number of contracts approaches infinity. Furthermore, our valuation 
operator is subadditive and satisfies a number of other appealing properties. Finally, our 
methodology allows us to decompose the value of any portfolio of pure endowment policies 
into a systematic component (due to uncertain aggregate mortality) and a non-systematic 
component (due to insuring a finite number of policies); see equation (4.55). 
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The remainder of this paper is organized as follows. In Section 2, we present our 
financial market, describe how to use the instantaneous Sharpe ratio to price the pure 
endowment, and derive the resulting partial differential equation that the price solves. 
In Section 3, we discuss qualitative properties of the risk-adjusted price from Section 2 
and show that it shares many properties with the (static) standard deviation premium 
principle (Gerber, 1979). In Section 4, we study properties of the price for n conditionally 
independent and identically distributed pure endowment risks. In particular, in Theorem 
4.11, we show that the price is subadditive with respect to n, and in Theorem 4.13, we 
show that the risk charge per person decreases as n increases. We also prove that if the 
hazard rate is deterministic, then the risk charge per person goes to zero as n goes to 
infinity (Theorem 4.20 and Corollary 4.21). Moreover, we prove that if the hazard rate is 
stochastic, then the risk charge person is positive as n goes to infinity, which reflects the 
fact that the mortality risk is not diversifiable (Theorem 4.20 and Corollary 4.22). Section 
5 concludes the paper. 

2. Instantaneous Sharpe Ratio 

In this section, we describe a pure endowment contract and present the financial 
market in which the issuer of this contract invests. We obtain the hedging strategy for the 
issuer of the pure endowment. We describe how to use the instantaneous Sharpe ratio to 
price the pure endowment and derive the resulting partial differential equation that the 
price solves. In Sections 3 and 4, we determine properties of the price. 

2.1. Mortality Model and Financial Market 

We begin with a stochastic model for mortality. We assume that the hazard rate A 
(or force of mortality) of an individual follows a diffusion process such that if the process 
begins at Ao > A for some positive constant A, then X t > A for all t e [0, T]. Thus, we 
require that the volatility of A goes to zero as A — > A from the right, and we require that 
the drift of A is positive for A close to A. The reason for requiring A to have a positive 
lower bound A will be apparent later. From a modeling standpoint, A could represent the 
hazard rate remaining (say, from accidents) after all biological causes of death have been 
removed. 

Specifically, we assume that 

d\ t = a(X t ,t)dt + b(t)(X t - X)dW t x , (2.1) 

in which W x is a standard Brownian motion on a probability space (O, JF, P). The volatility 
b is either identically zero, or it is a continuous function of time t bounded below by a 
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positive constant k on [0, T]. The drift a is a Holder continuous function of A and t for 
which there exists e > such that if < A — A < e, then a(A, t) > for all t £ [0, T]. 
After Lemma 3.3 below, we add additional requirements for a. Note that if b = 0, then A 
is deterministic, and in this case, we write X(t) to denote the deterministic hazard rate at 
time t. 

Suppose an insurer issues a pure endowment to an individual that pays 1 at time 
T if the individual is alive at that time. In Section 2.2, to determine the value of the 
pure endowment, we will create a portfolio composed of the obligation to pay this pure 
endowment and of default-free zero-coupon bonds that pay 1 at time T regardless of the 
state of the individual. Therefore, we require a model for bond prices, and we use a model 
based on the short rate and the bond market's price of risk. 

The dynamics of the short rate r, which is the rate at which the money market 
increases, are given by 



dr t = fJ,(r t , t)dt + a(r u t)dW t , 



(2.2) 



in which \i and a > are deterministic functions of the short rate and time, and W is a 
standard Brownian motion with respect to the probability space (O, JF, P), independent of 
W x . We assume that fx and a are such that r > almost surely. 

From the principle of no-arbitrage in the bond market, there is a market price of risk q 
for the bond that is adapted to the filtration generated by W; see, for example, Lamberton 
and Lapeyre (1996) or Bjork (1998). Moreover, the bond market's price of risk at time t 
is a deterministic function of the short rate and of time, that is, qt = q(rt,t). Thus, the 
time t price of a T-bond is given by 



F(r,t;T) = E Q 



- f T r s ds 




e J * 


r t =r 



(2.3) 



in which Q is the probability measure with Radon-Nikodym derivative with respect to P 
given by 



- J Q T q(r s ,s)dW s -± j^q 2 (r s ,s)ds 

dV 



(2.4) 



It follows that W Q , with W? = W t + f* q{r s ,s)ds, is a standard Brownian motion with 
respect to Q. 

From Bjork (1998), we know that the bond price F solves the following partial differ- 
ential equation (pde): 
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(2.5) 



in which \fi = fx — qa. In this paper, the horizon T is fixed, so henceforth we drop T from 
the notation of F. We can use this pde to obtain the dynamics of the bond price F(r s , s), 
in which we think of r t = r as given and t < s < T. Indeed, 

( dF(r s , s) = (r s F(r s , s) + q(r s , s)a(r s , s)F r (r s , s))ds + a(r s , s)F r (r s , s)dW s , 



As an aside, we could use other models commonly used in the literature and obtain 
the same conclusion that we reach after equation (2.17), namely, that we can factor the 
T-bond price from the "mortality price." We use a model involving the short rate and the 
bond market's price of risk for ease of presentation. 

2.2. Pricing via the Instantaneous Sharpe Ratio 

The insurer faces the unhedgeable risk that the individual's living or dying will be 
different from expected; therefore, the insurer demands a return greater than the sum of 
the return r on the money market and the "return" A from the mortality component. One 
measure of the risk that the insurer takes is the standard deviation of the change in the 
portfolio. A natural tie between the excess return and the standard deviation is the ratio 
of the former to the latter, the so-called instantaneous Sharpe ratio. In what follows, we 
find the hedging strategy to minimize the local variance of the change in the portfolio, then 
we set the price of the pure endowment so that the resulting instantaneous Sharpe ratio 
equals a given constant. We could set the instantaneous Sharpe ratio equal to a function 
of A and t, but we choose a constant for simplicity. 

The market for insurance is incomplete; therefore, there is no unique pricing mech- 
anism. To value contracts in this market, one must assume something about how risk is 
"priced." For example, one could use the principle of equivalent utility (see Zariphopoulou 
(2001) for a review) or the Esscher transform (Gerber and Shiu, 1994) to price the risk. 
We employ the instantaneous Sharpe ratio because of its analogy with the bond market's 
price of risk and because of the desirable properties of the resulting price. Because of these 
properties, we anticipate that our pricing methodology will prove useful in pricing risks in 
other incomplete markets. 

Denote the value (price) of the pure endowment by P = P(r, A, t), in which we explic- 
itly recognize that the price of the pure endowment will depend on the short rate r and 




(2.6) 
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the hazard rate A at time t. Suppose the insurer creates a portfolio IT with value U t at 
time t. The portfolio contains the obligation to pay the pure endowment at time T if the 
individual is alive at that time, namely —P. Additionally, the insurer holds 7r t T-bonds. 
Thus, II t = -P(r t , A t , t) + n t F(r t , t). 

By Ito's Lemma (Protter, 1995), the value of the portfolio at time t + h with h > 0, 
namely ILt+h, equals 

pt+h pt + h 

U t+h = U t - V^P(r s ,X s ,s)ds + a(r s ,s)(7i s F r (r s ,s)- P r (r s ,\ s ,s))dW s 



<t Jt 

-t+h pt+h 



pz-tn pT-tn 

- / b(s)(X s -X)P x (r Sl X Sl s)dW s x + / P(r s ,X s ,s)(dN s -X s ds) 
Jt Jt 
pt+h 

+ / n s (r s F(r s ,s) +q(r s ,s)a(r s ,s)F r (r s ,s))ds, 

(2.7) 

in which T> m , with m = m(r,t) a deterministic function of the short rate and time, is an 
operator defined on the set of appropriately differentiable functions on R + x (A, oo) x [0, T] 
by 

1 1 

V m v = v t + mv r + -a 2 v rr + av\ + -b 2 (X - X) 2 v X \ - Xv. (2.8) 

Also, in (2.7), N denotes a Poisson process with stochastic parameter A. Thus, IT jumps in 
value by P when an individual dies because the insurer is no longer responsible for paying 
1 at time T. 

In this single-life case, the process IT is "killed" when the individual dies. If we were 
to consider the price P^ for n conditionally independent and identically distributed lives 
(conditionally independent given the hazard rate), then would be a Poisson process with 
stochastic parameter nX such that IT jumps by P^ — p( n-1 ) when an individual dies. We 
consider P^ later and continue with the single-life case now. 

We next calculate the expectation and variance of Ht+h conditional on the information 
available at time t, namely Tt- First, 

pt+h 

E(Il t+h \F t ) = n - E r ' A -* / V^P(r s , X s , s)ds 

Jtth (2-9) 
+ E r ' A,t n s (r s F(r s , s) + q(r 8 , s)a(r s , s)F r (r s , s)) ds. 

Here, II t = IT is known at time t, and E r ' A,t denotes the conditional expectation given 
r t = r and At = A. Define the stochastic process Yh for h > by 
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pt+h pt+h 

Y h = U- V»P(r s ,X s ,s)ds + w s (r s F(r s ,s)+q( r s ,s)a(r s ,s)F r (r s ,s))ds. (2.10) 
Jt Jt 

Thus, E(Il t+h \F t ) = W^Yu, and from (2.7), we have 



cr(r s , s)(n s F r (r s , s) - P r (r s , X s , s)) dW s 

/t+h rt+h 
b(s)(X s - A)P A (r s , A s , s) dW s x + J P(r s , A s , s)(dN s - A s ds) 



H 

It follows that 



Var(n t+h |^ t ) = B((U t+h - EY h ) 2 \F t ) 



a 2 (r s , s)(n s F r (r s , s) - P r (r s , A s , s)) 2 ds (2J2) 

rt+h pt+h 

+ E r ^ / b 2 (s)(X s -X) 2 Pi(r s ,X s ,s)ds + E r ^ t / A s P 2 (r s , A s , s)ds. 
Jt Jt 



We choose nt in order to minimize the local variance lim/^o ^Var(ilt_|_^|jF t ), a dy- 
namic measure of risk of the portfolio; therefore, Tit = Pr{fti Xt,t)/F r (rt,t). Under this 
assignment, the drift and local variance become, respectively, 

lim i(E(n m |^) -H) = -V» Q P(r,X,t) + rP r (r,X,t)pp^, (2.13) 
h—>0 a ^(r, t) 

and 

lim ivar(n t+h |^ t ) = b 2 (t)(X - A) 2 P 2 (r, A, t) + AP 2 (r, A, t). (2.14) 

h— >o n. 

Now, we come to pricing via the instantaneous Sharpe ratio. Because the minimum 
local variance in (2.14) is positive, the insurer is unable to completely hedge the risk of the 
pure endowment contract. Therefore, the price should reimburse the insurer for its risk, 
say, by a constant multiple a of the local standard deviation of the portfolio. It is this a 
that is the instantaneous Sharpe ratio. 

From (2.14), we learn that the local standard deviation of the portfolio equals 



lim yivar(II t+/ ^) = y^ 2 (t)(A - A) 2 P 2 (r, A, t) + AP 2 (r, A, t). (2.15) 
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To determine the value (price) P, we set the drift of the portfolio equal to the short rate 
times the portfolio plus a times the local standard deviation. Thus, from (2.13) and (2.15), 
we have that P solves the equation 




r 



(2.16) 



for some < a < ^/X. Recall that II = -P + nF = -P + P r F/F r . It follows that 
P = P(r, A, t) solves the non-linear pde given by 



If we had been able to choose the investment strategy n so that the local standard 
deviation in (2.15) were identically zero (that is, if the risk were hedgeable), then the 
right-hand side of the pde in (2.17) would be zero, and we would have a linear differential 
equation of the Black-Scholes type. One can think of the right-hand side as adding a 
margin to the return of the portfolio because the pure endowment risk is not completely 
hedgeable due to the mortality risk. In addition to the unhedgeable mortality risk, the 
insurer also faces the somewhat diversifiable finite portfolio risk, that is, the risk that even 
if the hazard rate is deterministic, the actual number who survive until time T is different 
from expected. This risk is clearly present when selling a pure endowment to a single 
individual, but it is also present to some extent in any portfolio of finite size. In Section 
4.4, we decompose the risk loading in the price due to the finite portfolio risk and due to 
the stochastic mortality risk. 

If there were no risk loading, that is, if a = 0, then the price is such that the expected 
return on the price is r + A. The rate r arises from the riskless money market, and A arises 
from the expected release of reserves as individuals die. If a > 0, then the expected return 
on the price is greater than r + A. Therefore, a, the Sharpe ratio, measures the degree 
to which the insurer's total expected return is in excess of r + A, as a proportion of the 
standard deviation of the return. 

Before moving on to the following sections where we study properties of the solution 
of (2.17), we show that we can simplify P greatly. Indeed, P(r, A, t) = F(r, t)(p(\, i), in 
which F is the price of the T-bond and solves (2.5), and (p solves the non-linear pde 



P t + fi Q P r + \o 2 P rr + aP x + V(A - A) 2 P AA - (r + X)P 



= -a^b 2 (X-X) 2 P 2 + XP 2 
,P(r,A,T) = l. 



(2.17) 
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<p t + a^ x + l -b 2 {\ - A) Vaa -\<p= -a^62(A-A)Vi + A^ 2 , ^ 
¥>(A,T) = 1. 

The existence of a solution to (2.18) follows from standard techniques; see, for example, 
Walter (1970, Chapter IV, Section 36). A comparison result (see Section 3 of this paper) 
demonstrates that the solution is unique. 

The factorization P = Ftp is reminiscent of the standard actuarial method of pricing 
pure endowments in that ip represents the probability of paying the mortality-contingent 
claim, that is, the probability that the individual survives. We will see in Section 3.1 that 
we can interpret <p as a risk-adjusted survival probability. As mentioned at the end of 
Section 2.1, this factorization arises under other commonly used models for bond prices, 
as long as the risk driving the bond price (W in our case) is independent of the risk driving 
the stochastic hazard rate, W x . 

Consider the special case for which 6 = 0, that is, A is deterministic. Suppose X(t) is 
the solution of dX = a(A, s)ds with initial value Ao = A; then, (2.18) becomes the linear 
ordinary differential equation 

<p'(t) - (X(t) - ay/Mfj)<p(t) = 0, <p(T) = 1, (2.19) 

whose solution is 

<p(t) = e~ JT (A(S) - Q v 7 ^) ) ds ; (2.20) 

a type of probability of survival because we can think of X(t) - a^Xjt) > as a modified 
hazard rate. Indeed, note that if a = 0, then (2.20) is the physical probability that a 
person alive at time t survives to time T, and as a increases, (p increases. Therefore, we 
interpret (2.20) as a risk-adjusted probability of survival, in which a controls the degree 
to which we adjust (that is, increase) the physical probability of survival. When 6 ^ 0, a 
similar phenomenon occurs (see Theorem 3.8 below), and for this reason, we refer to the 
solution P of (2.17) as the risk-adjusted price for the pure endowment. 

Recall that we assume that < a < y/X. For the case of deterministic hazard, 
this implies that A(£) - a^/X{t) > for all t e [0,T], so that < <p(t) < 1 in (2.20). 
For stochastic hazard, we show in the next section that < <p(X,t) < 1 for all (A, t) G 
(A, oo) x [0, T]. It follows that, in general, F is an upper bound for the price. Observe 
that F is a natural upper bound for the price because it is the price we would charge if we 
knew the person could not die before T. 
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3. Qualitative Properties of the Risk- Adjusted Price 

In this section, we discuss qualitative properties of the risk-adjusted price P in (2.17) 
and show that it shares many properties with the (static) standard deviation premium 
principle (Gerber, 1979). To begin, we have the following proposition. 

Proposition 3.1. Suppose P c is the price, as determined by the method in Section 2, for 
a pure endowment with payment c > at time T if the individual is alive. Then, P c = cP, 
in which P is the risk-adjusted price for a payment of 1 at time T if the individual is alive. 

Proof. In the derivation of P for (2.17), it is clear that if we derive the price P c , then 
(2.17) still applies for determining P c with the terminal condition P c (r, A, T) = c. We can, 
then, write P c = cF(p = cP. □ 

Proposition 3.1 parallels the following well known fact concerning the (static) standard 
deviation premium principle H: If we define the standard deviation premium principle, as 
applied to a random variable X, by 



then H(cX) = cH(X) for c > 0. 

In what follows, we show that < P < F and P\ < 0, and we examine how the price P 
responds to changes in the model parameters. To this end, we need a comparison principle 
(Walter, 1970, Section 28). We begin by stating a relevant one-sided Lipschitz condition 
along with growth conditions. We require that the function g = g{X,t,v,p) satisfies the 
following one-sided Lipschitz condition: For v > w, 



with growth conditions on c and d given by 

0<c(A,£) < K(l + (ln(A-A)) 2 ), andO<d(A,t) < K(X - A)(l + | ln(A - A)|), (3.3) 

for some constant K > 0, and for all (A,t) G (A, oo) x [0, T]. Throughout this paper, 
we rely on the following useful comparison principle, which we obtain from Walter (1970, 
Section 28). 

Theorem 3.2. Let G = (A, oo) x [0,T], and denote by Q the collection of functions on G 
that are twice- differentiable in their first variable and once- differ entiable in their second. 



H(X) — EX + oVVarX, 



(3.1) 



g(X,t,v,p) - g(X,t,w,q) < c(X,t)(v - w) + d(X,t)\p - q 



(3.2) 
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Define a differential operator C on Q by 

Cv = v t + h?{t){\ - \fv xx + g(X, t, v, v x ), (3.4) 

in which g satisfies (3.2) and (3.3). Suppose v, w G Q are such that there exists a constant 
K > with v < e K ( ln ( A ~^) 2 and w > - e K( - ln ( x ~^ 2 for large X and for X close to X. Then, 
if (a) Cv > Cw on G, and if (b) v(X,T) < w(X,T) for all X > X, then v < w on G. 

Proof. Transform the variables A and t in (3.4) to y = ln(A — A) and r = T — t, and write 
v(y,r) = f(A,t), etc. Under this transformation, (3.4) becomes 

Cv = -v T + ^b 2 (r)vyy + h{y, r, v, v y ), (3.5) 

in which h(y, r, v, p) = — \b {r)p+g(y, r, v,p), and v is a differential function on R x [0, T]. 
Note that <p\ = e~ y (p y , so p = e~ y p in going from g to g. The differential operator in (3.5) 
is of the form considered by Walter (1970, Section 28, pages 213-215); see that reference 
for the proof of our assertion. 

The remaining item to consider is the form of the growth conditions in the original 
variables A and t. From Walter (1970), we know that analog of (3.2) and (3.3) for h are 

h(y,r,v,p) -h(y,r,w,q) < c(y,r)(v-w) + d(y,r)\p - q\, (3.6) 

with 

< c(y, t) < K(l + y 2 ), and < d(y, r)<K(l + \y\). (3.7) 

Under the original variables, the right-hand side of (3.6) becomes c(A, t)(v — w) +d(X, t)\p — 
q\, in which c(A, t) = c(y, r) and <i(A, t) = d(y, r)e y because p = e y p. Therefore, d{y, r) < 
K(l + \y\) becomes d(X,t) < Ke y (l + \y\) = K(X - A)(l + | ln(A - A)|). □ 

As a lemma for results to follow, we show that the differential operator associated 
with our problem satisfies the hypotheses of Theorem 3.2. 

Lemma 3.3. If we define g by 

g(X, t, v,p) = a(X, t)p -Xv + a^b 2 (t)(X - A) V + Xv 2 , (3.8) 

then g satisfies the one-sided Lipschitz condition (3.2) on G. Furthermore, if \a(X, t)\ < 
K{X - A)(l + | ln(A - A)|), then (3.3) holds. 

Proof. Suppose v > w, then 



13 



g(X,t,v,p) - g(X,t,w,q) = a(X,t)(p- q) - X(v - w) 



+ a | y/b 2 (t)(X- A)V + Xv 2 - \Jb 2 (t)(X - X) 2 q 2 + Amj 2 j 

< (\a(X,t)\ + ab(t)(X-X))\p-q\ - (x - aVx) (v - w) 
<(\a(X,t)\ + ab(t)(X-X))\p-q\. 



(3.9) 



Recall that a < y/£. Also, we use the fact that if A > B, then VC 2 + A 2 - VC 2 + B 2 < 
A — B, as we demonstrate below in Lemma 4.5. Thus, (3.2) holds with c(A, t) = and 
d{\t) = \a(X,t)\ + ab(t)(X - A). Note that d satisfies (3.3) if |a(A, *)| < K(X - A)(l + 



Assumption 3.4. Henceforth, we assume that the drift a satisfies the growth condition 
in the hypothesis of Lemma 3.3. For later purposes (for example, see Theorem 3.7), 
we also assume that ci\ is Holder continuous and satisfies the growth condition \a\\ < 



In the next two subsections, we apply Theorem 3.2 and Lemma 3.3 repeatedly to 
determine qualitative properties of the risk-adjusted premium P. 

3.1. Interpreting <p as a Survival Probability 

In our first applications of Theorem 3.2, we show that ip shares two important prop- 
erties with physical survival probabilities, namely that < f < 1 and ip\ < 0. In fact, we 
have an even tighter upper bound on <p by observing that A is a lower bound for A. 

Theorem 3.5. < ip(X,t) < e"^-^)^-*) f or (A,£) e G. 

Proof. Define the differential operator C on Q by (3.4) with g given in (3.8). Because ip 
solves (2.18), we have Cp> = 0. Also, 



Because £ e -(A-«VA)( T -*) < Cip and <p(\,T) = 1 = e -^- a ^)( T - T \ Theorem 3.2 and 
Lemma 3.3 imply that v? < e "( A " Q v / A)( T - t ) on G. 

Similarly, denote by the function that is identically on G; then, £0 = = Lp>. 
Because additionally (p(X, T) = 1, Theorem 3.2 and Lemma 3.3 imply that < <p on G. □ 



ln(A-A)|). 



□ 



K(l + (ln(A-A)) 2 ). 




(3.10) 
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The upper bound in Theorem 3.5 is tight. Indeed, suppose A = A + e for some constant 
e > 0. Then, the solution to (2.18) is given by e "( ( - +e) " Q v / ^)( T - t ), which can be made 
arbitrarily close to e~^-~ a ^^ T ~ t ^ > (uniformly on [0, T]) by choosing e small enough. 

We have the following corollary of Theorem 3.5 that gives us natural bounds for the 
risk-adjusted price. F is a natural bound for the price because it is the price of a default- 
free bond, that is, a bond that pays regardless of whether the individual is alive. 

Corollary 3.6. < P(r, A, t) < F(r, t) for (r, A, t) G R + x G. 

Proof. Because P(r, A, t) = F(r, t)ip(X, t) and e~^-~ a ^ {T ~ t] < 1, the result is immedi- 
ate from Theorem 3.5. □ 

We end this subsection with a proof that (p\ < 0. This result is intuitive for physical 
survival probabilities in that if the current hazard rate A increases, then the probability of 
surviving until time T decreases. 

Theorem 3.7. (p\(X,t) < for (\,t) G G. 

Proof. To prove this assertion, we apply a modified version of Theorem 3.2 to the special 
case of comparing (p\ with the zero function 0. From Walter (1970, Section 28, pages 
213-215), we see that we only need to verify that (3.2) holds for v > = w = q. First, 
differentiate y's equation with respect to A to get an equation for f = 

' ft + (ax-X)f+(a + b 2 (X - A))/ A + \b 2 {\ - A) 2 /aa - V 

= a fr 2 (A - A)/ 2 + 6 2 (A - A) 2 /A + y 2 + \<pf (3 n) 

V& 2 (A-A) 2 / 2 + Av? 2 

(f(X,T) = 0. 

Define a differential operator C on Q by (3.4) with g given by 

2 , 6 2 (A-A)z; 2 + 6 2 (A-A) 2 ^+i^ 2 + A^ 



g(\,t,v,p) = (a\- X)v + (a + b (A- X))p — (p + a 



^b 2 {X- X) 2 v 2 + X^ 2 

(3.12) 

To apply a modified version of Theorem 3.2, verify that (3.2) and (3.3) hold for v > = 
w = q. It is not difficult to show that in this case, 

g(X,t,v,p) - g(X,t, 0,0) < (\a x \ + ab - (A - aVX)v + (\a\ + b 2 (X - A) + ab(X - X))\p\ 

< (\a x \ + ab)v + (\a\ + b 2 (X - A) + ab(X - X))\p\. 

(3.13) 
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Thus, by Assumption 3.4, g satisfies (3.2) with the corresponding c and d satisfying the 
growth conditions in (3.3). 

Next, note that because f = <p\ satisfies (3.11), Cf = 0. Also, CO = (p(—l + 
a/{2\/X)) < because A > A > a 2 . These observations, together with /(A, T) = 0, 
imply that / = <p\ < on G. □ 

As in the relationship between Theorem 3.5 and Corollary 3.6, we have the immediate 
corollary of Theorem 3.7 that P\ < on G. 

3.2. Comparative Statics for P 

Our next results show that as we vary the model parameters, the price P responds 
consistently with what we expect. 

Theorem 3.8. Suppose < a.\ < a 2 < y/X, and let P ai be the solution to (2.17) with 
a = a*, fori = 1,2. Then, P Ql (r, X,t) < P a2 (r,X,t) for all (r, A, t) G R + x G. 

Proof. Because F in (2.5) is independent of a, it is enough to show that (p ai < (p a2 on 
G, in which ip ai has the obvious meaning. Define a differential operator £ on Q by (3.4) 
and (3.8) with a = a±. Because </? ai solves (2.18) with a = «i, we have £</? ai = 0. Also, 



= + a<p? + -6 2 (A - A)Vaa - A^ 2 + a iy Jb?(\- m<P?) 2 + H<P°») 2 
= -(a 2 - ai)y / 6 2 (A-A) 2 (^ 2 ) 2 + A(^)2 < = Cl p a K 

(3.14) 

In addition, both (f ai and cp a2 satisfy the terminal condition (p ai (A, T) = 1. Thus, Theorem 
3.2 and Lemma 3.3 imply that </? ai < </?° 2 on G. □ 

Theorem 3.8 states that as the parameter a increases, the price P a increases; this 
result justifies the use of the phrase risk parameter when referring to a. It is clear that 
the price obtained from the standard deviation premium principle, as defined in (3.1), also 
increases with a. We have the following corollary to Theorem 3.8. 

Corollary 3.9. Let P a0 be the solution to (2.17) with a = 0; then, P a0 < P a for all < 
a < y/X, and we can express the lower bound P a0 as follows: P Q °(r, A, t) = F(r, % a0 (A, t), 
in which </? a0 is given by the physical probability of survival, namely 



V^°(A,t) = E 



e J t 



J \ s ds 



Xt = X 



(3.15) 



where X s follows the process given in (2.1). 
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Proof. Theorem 3.8 implies that P a0 < P a for all < a < y/X, and by substituting 
a = in (2.18), the Feynman-Kac Theorem (Karatzas and Shreve, 1991) implies the 
representation of (p a0 in (3.15). □ 

Corollary 3.9 justifies the use of the phrase risk-adjusted price when referring to P a 
because P a > P Q °, and we can interpret _P Q ° as a risk-neutral price due to the fact that P a0 
equals the product of the bond price and the probability of paying the pure endowment. 
We call P a — P a0 the risk charge that compensates the insurer for both the finite portfolio 
risk and the stochastic mortality risk. In Section 4.4, we decompose the risk charge into 
these two components after we study the price for a portfolio of n pure endowment risks. 

Next, we examine how the risk-adjusted price P varies with the drift and volatility of 
the stochastic hazard rate. 

Theorem 3.10. Suppose ai(X,t) < a 2 (A,t) on G, and let P ai denote the solution to (2.17) 
with a = at, fori = 1,2. Then, P ai (r, X,t) > P a2 (r, X,t) for all (r, A, t) G R + x G. 

Proof. Because F in (2.5) is independent of a, it is enough to show that f ai > f a2 , in 
which <p ai has the obvious meaning. Define a differential operator C on Q by (3.4) and 
(3.8) with a = a±. Because <p ai solves (2.18) with a = a\, we have £</? ai = 0. Also, because 
V?^ 2 < 0, then 



= + aw? + -b 2 (X - X)\Z - Xip a2 + a^JbHX-X)H^r + X(^r (g J6) 
= -(a 2 - ai )<p a x 2 >Q = &p a K 

In addition, both (p ai and p a2 satisfy the terminal condition </? ai (A,T) = 1. Theorem 3.2 
and Lemma 3.3 imply that tp ai > ip a2 on G. □ 

The result parallel to Theorem 3.10 for the standard deviation premium principle in 
(3.1) is that H decreases as EX decreases. 

Theorem 3.11. Suppose < bi(t) < 62^) on [0, T], and let P bi denote the solution to 
(2.17) with b = 6,, for i = 1, 2. // P* x > for i = 1 or 2, then P 6l (r, A, t) < P ba (r, A, t) 
for all (r, A, t) E R + x G. 

Proof. Because F in (2.5) is independent of 6, it is enough to show that (p bl < (p b2 if 
Vaa ^ 0> m which ip bi has the obvious meaning. Suppose ip b x \ > 0, and define a differential 
operator C on Q by (3.4) and (3.8) with b = b\. Because (p bl solves (2.18) with b = 61, we 
have &p bl = 0. Also, 
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V 2 = ^ + V A 2 + \bl(X - X)Vx\ - A/ 2 + a^bl(X-X)H^) 2 + H^) 2 
= -U b l~ b i) ( A _ A) Vaa 



(3.17) 



<0 = Cip bl . 

In addition, both (p bl and (p b2 satisfy the terminal condition (p bi (X,T) = 1. Thus, Theorem 
3.2 and Lemma 3.3 imply that (p bl < (p b2 on G. In a similar fashion, we obtain the same 
result if ip b x \ > 0. □ 

From Theorem 3.11, we see that if ip is convex with respect to A, then the risk-adjusted 
price increases if the volatility on the stochastic hazard rate increases. The parallel result 
for the standard deviation premium principle in (3.1) is that H increases as VarAT increases. 



4. The Risk-Adjusted Price for a Portfolio of Risks 

In this section, we study properties of the price P^ = P^ (r, A, t) for n condition- 
ally independent and identically distributed pure endowment risks. First, we present the 
equation that P^ solves. In Sections 4.1 and 4.2, we parallel the results of Sections 3.1 
and 3.2, respectively. In Section 4.3, we show that P^ is subadditive. In Section 4.4, 
we show that the risk charge per person decreases as the n increases. We also show that 
if the hazard rate is deterministic, then the risk charge per person goes to zero as n goes 
to infinity. Moreover, we show that if the hazard rate is stochastic, then the risk charge 
person is positive as n goes to infinity, which reflects the fact that the stochastic mortality 
risk is not diversifiable. 

As discussed in the paragraph preceding equation (2.9), when an individual dies, the 
portfolio value IT jumps by P^ — p( n_1 ). By paralleling the derivation of (2.17), one can 
show that solves the non-linear pde given by 

' P t {n) + Li Q P^ n) + \y 2 P^ + aP ( x n) + \b 2 (X - A) 2 Pff - rPW - nX (p^ - P^ 1 )) 

= -a^b 2 (X-X) 2 (p| n) ) 2 + nA (p(-) _ p(--i)) 2 
P^(r,X,T) = n. 

(4.1) 
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The initial value in this recursion is P(°) = 0, and the price P as defined by (2.17) is pW. 

As in Section 2, we can multiplicatively separate the variables r and A in p( n \ Indeed, 
p( n )(r, X,t) = F(r, t)ip { - n \X, t), in which F solves (2.5) and ip^ solves the recursion 

' + + \b\X - A) VS " n\ (y™ - <p^) 

= -a^b 2 (X - A) 2 (<^i n) ) 2 + nX (pW - v^" 1 )) 2 , (42) 

y*>(A,T)=n, 

with initial value = 0. Note that </? in (2.18) equals (p^ 1 '. 

Throughout this section, we apply Theorem 3.2 with g = g n defined by (4.3) below. 
We have the following lemma whose proof we omit because it parallels that of Lemma 3.3. 

Lemma 4.1. Define g n , for n > 1, by 

g n (X, t, v,p) = a(A, t)p - nX (v - v? (n_1) ) + a\]b 2 {X - A) V + nX (v - v?^" 1 )) 2 , (4.3) 

in which ip( n ~ 1 ' ) solves (4.2) with n replaced by n — 1. Then, g n satisfies the one-sided 
Lipschitz condition (3.2) onG. Furthermore, if\a(X, t)\ < K(X — A)(l + | ln(A — A)|), then 
(3.3) holds. 



4-1- Interpreting ^V 7 ^ as a Survival Probability 

In the first application of Theorem 3.2 and Lemma 4.1, we show that < ±v? (n) < 

e -(±-*VX)(T-t) forn > L 

Theorem 4.2. < ^ n) {X,t) < ne"^ -0 ^)^-*) for (X,t) G G and for n> 0. 

Proof. We proceed by induction to prove the upper bound. For the ease of presentation 
in this proof, define h(t) = e~(- _a v / ^)( T_ ^. It is clear that the inequality holds for n = 
because cp^ = 0. For n > 1, assume that < <^( n_1 ) < (n — l)h(t), and show that 
< ip {n) < nh(t). 

Define the differential operator C on Q by (3.4) with g = g n from (4.3). Because (p( n > 
solves (4.2), we have C<p^ = 0. Also, 

Cnh(t) = n(X - ay/X)h(t) - (nX - av^A) (nh(t) - t^™" 1 ^ 

< n(X - a^X)h(t) - (nX - aVnX)(n - (n - l))h(t) (4.4) 
n(X - a VA) - (nX - aVnX)] h(t) < 0. 
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Because ^ne"^" ^)^"*) < &p^ and ip^(X,T) = n, Theorem 3.2 and Lemma 4.1 
imply that v? (n) < ne~(- _a v^) (T_t) on G. 

Similarly, we prove the lower bound via induction. Suppose that </?( n_1 ) > for n > 1, 
and show that (p^ > 0. Let C be the differential operator from the first part of this proof, 
and denote by the function that is identically on G; then, CO = (nX+aVnX)^ 71-1 ^ > 0. 
Because CO > C<pW and ^ (n) (X, T) = n, Theorem 3.2 and Lemma 4.1 imply that < ip^ 
on G. □ 

An immediate corollary of Theorem 4.2 is that < P^ n ' < nF, as in Corollary 3.6, 
so, once again the price per risk ^P^ lies between and F. We end this subsection by 
extending Theorem 3.7, namely ip\ < 0, to an arbitrary number of risks. We first present 
a lemma which we use in its proof. 

Lemma 4.3. (A, t) > cp^'V (A, t) for (A, t) e G and for n>l. 

Proof. We proceed by induction. This inequality is true for n = 1 because ip^ = ip > 
= <p(°) by Theorem 3.5. For n > 2, assume that </?( n_1 ) > ip( n ~ 2 \ and show that 

^(n) > y,(n-l)_ 

Define a differential operator £ on £ by (3.4) with g = g n from (4.3). Because <^ n ) 
solves (4.2), we have Cip^ = 0. Also, 

£^(n-D = („ _ 1)A ^(n-1) _ ^(n-2) j _ nX ^(n-1) _ ^(n-1) j 

+ a^6 2 (A - A) 2 (v?i n_1) ) 2 + nA (v?^" 1 ) - v^" 1 )) 2 



(4.5) 



- a\J b 2 (X - A) 2 (v?i n_1) ) 2 + (n-l)A(^(-" 1 ) -^(--s)) 2 
> (n - 1)A (V^ - ^ n ~V) + ab(X - A) U n_1) 

- a (b(X - A) |^ n - 1} | + y/(n-l)X (V^ - ^ n ~ 2) ) ) 

= ((ra - 1)A - ay/(n - 1)a) (^ (n_1) - v? (n ~ 2) ) > = 

Note that the first inequality follows from the subadditivity of the square root, that is, 
vM 2 + B 2 < \A\ + \B\. We also use the induction hypothesis that ip^ n ^ > ip^ n 2 \ In 
addition, <p (n) (^ T) = n> n-1 = v? (n-1) (-^ T). Thus, Theorem 3.2 and Lemma 4.1 imply 
that v? (n) > on G. □ 

Lemma 4.3 is interesting in its own right because it confirms our intuition that P^ n ' 
increases with the number of policyholders n. 

Theorem 4.4. ip^ (A, t) < for (A, t) e G and for n>0. 
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Proof. We proceed by induction. We know that the inequality holds when n = and 1. 
For n > 2, assume that </^ n_1 ' ) < 0, and show that tp^ < 0. As in the proof of Theorem 
3.7, to prove this assertion, we apply a modified version of Theorem 3.2 to the special 
case of comparing cp^ with the zero function 0. First, differentiate </?( n )'s equation with 
respect to A to get an equation for = <p x n ' 



' f^ + (a x - *A)/<»> + nXf^ + (a + b\\ - A))/<"> + '-b\\ - A) 2 /, 
-n(^)-^- 1 )) 

b 2 (X - A) (/^)) 2 + 6 2 (A - X) 2 f^f{ n) + \n (<p( n ) - D) 2 

y^A - A) 2 (/W) 2 + nX (^0 - v?^" 1 )) 2 
n A (^W - ^("-i)) (/(») - 

y^A - A) 2 (/W) 2 + nA (</» - ^(--i)) 2 ' 
k /< n >(A,T) = 0. 

Define a differential operator £ on Q by (3.4) with g = g n given by 



g n (X, t, v,p) = (a x - nX)v + nXf^ + (a + b 2 (X - X))p - n (v? (n) - 

+ a b2 ( X ~ A)^ 2 + b 2 (X- Xfvp + | (v? (n) - V7 (n_1) ) 2 + nX (ip^ - ip( n -V) (v - f n -V) 

^b 2 (X-X) 2 v 2 + nX (^0 -^(--i)) 2 

(4.7) 

From Walter (1970, Section 28, pages 213-215), we see that we only need to verify that 
(3.2) holds for v > = w = q. It is not difficult to show that Lemma 4.3 implies that 
g n satisfies the inequality in (3.13). Thus, by Assumption 3.4, g n satisfies (3.2) with the 
corresponding c and d satisfying the growth conditions in (3.3). 

Next, note that because f^ = ^>[ n) satisfies (4.5), Cf^ = 0. Also, we have CO = 
(nX — aVnXj /( n_1 ) — — a/2^Jn/ X^j [tp^ — </?( n_1 )) < by the induction assumption, 
by Lemma 4.3, and because A > A > a 2 . These observations, together with f^ n \X,T) = 0, 
imply that = p[ n) < on G. □ 

4-2. Comparative Statics for P( n ) 

In this section, we present properties of P^ to parallel those in Section 3.2. We show 
that as we vary the model parameters, the price P^ responds consistently with what we 
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expect, as did P = P^ 1 ' in Section 3.2. We begin with a lemma that will help in proving 
Theorem 4.6 below. 

Lemma 4.5. Suppose A> B andC are constants; then, \JC 2 + A 2 < (A—B)+\ / C 2 + B 2 . 

Proof. By squaring both sides, we see that this inequality is equivalent to C 2 +A 2 < (A 2 — 
2AB+B 2 ) + (C 2 +B 2 )+2(A-B)VC 2 + B 2 , which simplifies to < (A-B)(VC 2 + B 2 -B), 
which is clearly true. □ 

Theorem 4.6. Suppose < «i < a 2 < \fX, an d ^ P ( ^ n ' > ' ai be the solution to (4.1) 
with a = oti, for i = 1,2 and for n > 0. Then, P^' ai (r,X,t) < PW« (r, A, t) for all 
(r, X,t) G R + x G. 

Proof. It is enough to show that <^>( n )' ai < <^( n )> a 2 on G, in which <£>( n )> ai has the obvious 
meaning. We proceed by induction. It is clear that the inequality holds for n = because 
(p°' ai = for i = 1,2. For n > 1, assume that <^( n-1 )' ai < ^( n_1 )' a2 , and show that 

^(n),ai < ^(n),a 2 _ 

Define a differential operator £ on ^ by (3.4) with g = g n from (4.3) with a = cx\. 
Because <£>( n )> ai solves (4.2) with a = a±, we have C(p^ ,ai = 0. Also, 



£^(n),«2 = _ nA ^(n),aa _ ^(n-l),^ + nA ^(n),a 2 _ ^(n-l),a 2 ^ 

+ ai\Jb 2 (X - A) 2 (^ n) ' a2 ) 2 + nA (^(n),^ _ v? (n-l),a 1 ) 2 

- a 2 \jb 2 {X - A) 2 (^ n) ' a2 ) 2 + ?lA (^(n),^ _^(n-l),a 2 ) 2 

= -nA (V™- 1 )' 02 - </? (n_1) ' ai ) 
+ ai |^6 2 (A- A) 2 (^ n) ' a2 ) 2 + nA (v?H'« 2 - ^MWf 



(4.8) 



- y 6 2 (A - A) 2 (<^i n),a2 ) +nA(v?^' a2 -v? (n - 1} ' a2 )' 

- (a 2 - ai)^ 2 (A- A) 2 (<^ n) ' a2 ) 2 + nA - ^(n-i),^) 2 

< - (nX - aiVnA) (V™ -1 ^ 2 - v^" 1 ^ 1 ) < = £<^' a \ 

Note that the first inequality follows from Lemma 4.5 with A = y/nX (ip( n ^ a2 — <^( n ~ 1 )' Q!1 ), 
B = ^n~X (cp^^ 2 - t^" 1 )^ 2 ), and C = b{X - A)<^ n),a2 , from the induction hypothesis, 
and from a 2 > ot\. In addition, both ip( n ^> ai and ip( n ^ a2 satisfy the terminal condition 
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p(»).«i(A,T) = 1. Thus, Theorem 3.2 and Lemma 4.1 imply that ^ n )^ < ipM^ on G. 
□ 

Theorem 4.6 extends Theorem 3.8 and states that as the parameter a increases, the 
risk-adjusted price _p( n )' a increases. We have the following corollary to Theorem 4.6. 

Corollary 4.7. Let p( n )> a0 be the solution to (4.1) with a = 0; then, p( n )> a0 < p(»).« for 
all < a < y/X, and we can express the lower bound p(™)> a0 as follows: P^> a0 (r,\,t) = 
nF(r, t)</? a0 (A, t), in which </? a0 is given by (3.15). 

Proof. It is straightforward to show that mp a0 solves (4.2) with a = 0, and the result 
follows. □ 

Note that rup a0 is the expected number of survivors under the physical measure, so the 
lower bound of ^P^ (as a approaches zero) is the same as the lower bound of P, namely, 

Next, we examine how the risk-adjusted price p( n ) varies with the drift and volatility 
of the stochastic hazard rate. We state the following two theorems without proof because 
their proofs extend those of Theorems 3.10 and 3.11, respectively, as the proof of Theorem 
4.6 extends the one of Theorem 3.8. 

Theorem 4.8. Suppose ai(A, t) < a 2 (A, t) on G, and let p( n )> ai denote the solution to (4.1) 
with a = a u fori = 1,2. Then, P( n '' ai (r, A, t) > P^> a2 (r,\,t) for all (r, A, t) e R + x G. 

Theorem 4.9. Suppose < b\{t) < 62^) on [0, T], and let p( n )' &i denote the solution 
to (4.1) with b = b t , for i = 1,2. // > for i = 1 or 2, then P^' bl (r, A,t) < 

p(n),b 2 ^ A ^ Q j Qr aU ^ A ^ G R + x G. 



4.3. Subadditivity of 

Subadditivity holds for the standard deviation premium principle, that is, H(X) + 
H(Y) > H(X + Y); thus, we expect it to hold for our pricing rule. We next show that 
P^ is subadditive. Specifically, we show that for m, n nonnegative integers, the following 
inequality holds: 

p(m) _^ p(n) > p(m+n)^ ^ g-j 

Subadditivity is a reasonable property because if it did not hold, then buyers of insurance 
could insure risks separately and thereby save money. 

We begin with a lemma that we will help us prove (4.9). 
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Lemma 4.10. Suppose A > C > B , B\, and C\ are constants; then, for nonnegative 
integers m and n, 

yJ{B x + C x ) 2 + (to + n)A 2 - ^{A - C) < ^ B 2 + mB 2 + ^C 2 + nC 2 + V^(A - B). 

(4.10) 

Proof. The left-hand side of (4.10) is nonnegative because A > A — C > 0. Square both 
sides of (4.10) to get 



(B x + C x ) 2 + (to + n)A 2 + n(A - C) 2 - 2^(A - C) ^{B x + C x ) 2 + (m + n)A 2 
<B\ + mB 2 +C 2 X + nC 2 + m(A - B) 2 + 2^(A - B) j ^B 2 + mB 2 + ^C 2 + nC 2 | 

+ 2sJb 2 + mB 2 ^C 2 + nC 2 , 

(4.11) 

which simplifies to 



B x Cx + nA(A - C) + mB(A - B) < ^{A - C)^(B X + C x ) 2 + (m + n)A 2 
+ Vn~i(A - B) ^Jb 2 + mB 2 + y ' C 2 + nC 2 | + ^ B 2 + mB 2 ^C 2 + nC 2 . 

If the left-hand side of (4.12) is nonpositive, then we are done. Suppose the left-hand 
side of (4.12) is positive, so that after squaring both sides, this inequality is equivalent to 

B\C\ + 2B x Cx {nA(A - C) + mB(A - B)} + n 2 A 2 (A - C) 2 

+ 2mnAB(A -B)(A-C)+ m 2 B 2 (A - B) 2 
< n(A - C) 2 {(B x + Cx) 2 + (m + n)A 2 } + m(A - B) 2 {B 2 +C 2 + mB 2 + nC 2 } 
+ (B 2 +mB 2 )(C 2 + nC 2 ) + D, 

(4.13) 

in which D > is a sum of nonnegative square- root terms. Inequality (4.13) simplifies to 

< to {C 2 (A - B) 2 + (B X (A - B) - C X B) 2 } + n {B 2 (A - C) 2 + (B X C - C X (A - C)) 2 } 
+ D + mn{(A(A - C) - C(A - B)) 2 + B 2 C 2 + 2A(C - B)(A- B)(A- C)} , 

(4.14) 

which is true because C > B. □ 
Theorem 4.11. If m and n are nonnegative integers, then p( m ) + p( n ) > p( m + n ) . 
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Proof. We prove this inequality by induction on the sum m + n. We know that (4.9) 
holds when m or n equals or when m + n equals or 1 because P^ = 0. For m + n > 2 
with m > 1 and n > 1, assume that + > P (k+e '> for all nonnegative integers k 
and £ such that k + £<m + n-l. We proceed to show that P^ + P^ > P (m+n) . 
The function £ = + (p( n ^ solves 



' £t + <A + ^ 2 (A - A) 2 6a - mA (V m) - cp^- 1 ^ - nX (tp™ - ^ n ~^ 
= -a\Jb 2 (X-X) 2 (^ m) ) 2 +mA (^M - ^-i)) 2 

- a^b 2 (X - A) 2 (vi n) ) 2 + nX - y^" 1 )) 2 , 

U(A,T)=m + n, 

and = <^( m+n ) solves 

' & + a0 A + l -b 2 {X - A) 2 AA - (m + n)A (0 - y^"" 1 )) 

= -a\Jb 2 (X - A) 2 2 + (m + n)A (0 - ^(-+--i)) 2 , ( 4 - 16 ) 
k 0(A,T) = m + n. 

Define a differential operator £ on ^ by (3.4) with g = g m +n from (4.3). From (4.16), 
we deduce that Ccfi = 0. Also, because £ solves (4.15), we have 

Ci = mX (V m) - ^ (m_1) ) + nX (V n) - ^ (n - 1} ) - (m + n)X - ^ m + n ~ 1 ^ 



b 2 (X 


-A) 2 (^i m) ; 


2 2 

+ mA (v?( m ) - ^( m ~ 1 )) 


b 2 (X 


-A) a (^ B) ) 


2 2 

+ nA (y?( n ) - ^("-i)) 



(4.17) 

+ a\Jb 2 (X - A) 2 £ 2 + (m + n)A (f - ^(m+n-i))^ 

< - (V m) + ip (n ~ 1} - ^ (m+n - 1} ) (nA - avmA) 

< = C<j). 

The last inequality follows from the induction hypothesis and from A — aV~X > for 
A > A > a 2 . The first inequality in (4.17) follows from Lemma 4.10 with the assignments 
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B x = b(X- A)^ m) , C x = b(\-\)<p[ n \ A = V\ (£ - p("H-»-i)), S = - t^™" 1 )), 

and C = \f\ (ip( n ^ — </?( n_1 )). By the induction assumption, we have A > B and A > C; 
without loss of generality, C > B. The functions £ and <fi satisfy the same terminal 
condition when t = T. Thus, Theorem 3.2 and Lemma 4.1 imply that £ > 0, or equivalently 

p(m) _|_ p(n) > p(m+n)_ |-| 

Theorem 4.11 gives us another proof - alternative to the one in Theorem 4.2 - that 
p(n) ^ namely, < nP^ < nF : in which we use Corollary 3.6 to assert that 

pi 1 ) = P < F. 

4-4- Limiting Behavior of P^ 

We next consider the limiting behavior of P( n ) . To motivate the results of this section, 
reconsider the standard deviation premium principle from (3.1) as applied to pricing n risks 
Xi, X 2 , • • • , X n that are identically distributed to a random variable X and conditionally 
independent given the random hazard rate. For concreteness, suppose that Xi is the 
indicator random variable of the event that individual i will be alive at time T. Then, 
Xs = Yh=i Xi is the total number of survivors at time T, so that 'E(Xs) = nEX, and 

Var(X 5 ) = Var[E(X s |5)] + E[Var(X s |5)] 

(4.18) 

= n 2 Var[E(X\S)] + nE[Var(X\S)], 

in which S is the cr-algebra generated by W x , the Brownian motion driving the stochastic 
hazard rate. Therefore, the static standard deviation premium principle (3.1) gives us 



H{X S ) = riEX + ctn\IVar[E(X\S)] + -E[Var(X\S)]. (4.19) 

V n 

Note that ^H(Xs) decreases as n increases. Also, in the limit, we have 

lim -H(Xs) = EX + aJ\ar[E(X\S)], (4.20) 

n^oo ft 

which is strictly greater than EX if a > and if Var[E(X|<S)] > 0. This last inequality will 
hold if b is uniformly bounded below by k > 0. Otherwise, if b = 0, then Var[E(X|«S)] = 0, 
and the right-hand side of (4.20) is simply EX. 

In this section, we show that P( n ) behaves much as H(Xs) does in (4.19) and (4.20). 
In Theorem 4.13, we show that the price per risk, ^pi n \ decreases as n increases; that 
is, by increasing the number of individuals insured, we reduce the risk per individual (as 
measured by the price). This result is consistent with what we expect, as inspired by 
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(4.19). The question answered by the results that follow Theorem 4.13 is how far does 
decrease, and we obtain results in close parallel to (4.20). 
We begin with a useful lemma. 

Lemma 4.12. If n > 2, and if A > C > and B\ are constants, then the following 
inequality holds 



B\ + -C 2 < Vn~^2(A — C) + \ l B\ + —^—{{n - 1)C - (n - 2)A) 2 . (4.21) 
n V n — 1 

Proof. By squaring both sides of (4.21), we can show that (4.21) is equivalent to 

-C 2 <(n-2)(A- Cf + —^—((n - 1)C - (n - 2)A) 2 

n n — 1 . t . 

— (4.22) 



+ 2\fn~^2{A - C)^B~l + ~ 1)C - (n - 



2)Ay 



Because B\ is arbitrary, (4.22) is true if and only if (4.22) holds when B\ = 0, which is 
equivalent to the following after taking the square root of the resulting right-hand side of 
(4.22): 



n - 2 \ , , _ „ f , 7i-2 1 



Inequality (4.23) holds if y/n - 1 - - ^ > 0, which is true for n > 2. □ 

Theorem 4.13. ^P^ decreases with respect to n > 1. 
Proof Define C (n) = ^ {n) for n > 1. Note that C (n) solves 

' Ci n) + aCi n) + \b\\ - X) 2 CS - A (nC (n) - (n - l)^) 



= -axJb 2 (X - X) 2 (Cf)' + «W - (n - l)^- 1 )) 2 , (424) 
[C (n) (A,T) = l, 
with C (1) = <£> (1) = (p. 

We proceed by induction and show that < C( n_1 ) for n > 2. We first show that 
£( 2 ) < q 1 ) = p. Define a differential operator C on Q by (3.4) with g = g 2 given by 



g 2 {X, t, v,p) = a(A, t)p - X(2v - <p) + a\/b 2 (X - A) V + ^X(2v - v?) 2 . (4.25) 
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Clearly, g 2 satisfies (3.2) and (3.3), so we can apply Theorem 3.2. Note that because 
solves (4.24) with n = 2, £C (2) = 0. Also, 



C<p = a |^ 2 (A - A) VI + ~ \Jb 2 (\ - A) VI + A^ 2 | < = £( {2 \ (4.26) 

Additionally, C (2) (A,T) = 1 = ip(X,T); thus, C (2) < C (1) = <P- 

Next, assume that for n > 3, C (n_1) < C (n ~ 2 \ and show that C (n) < C (n_1) - Define a 
differential operator V on Q by (3.4) with g = g n given by 



— a 



£ n (A,t, = ap-A (m; - (n - lR (n_1) ) + a^jb 2 (X - X) 2 p 2 + ^A (raw - (n - l)^"- 1 )) 2 . 

(4.27) 

Clearly, g n satisfies (3.2) and (3.3), so we can apply Theorem 3.2. Note that because 
solves (4.24), £>C (n) = 0. Also, 

= (n- 2)A (C (n_1) - C (n " 2) ) 

jyW - X) 2 (Ci n - 1} ) 2 + ^A ((n - fR^" 1 ) - (n - 2R(- 2 )) 2 

-^^-^(cr^v^A^- 1 )) 2 ! 

< ((n - 2)A - a v / (^-2)A) (c (n_1) - C (n " 2) ) < = ©C (n) - 

(4.28) 

The first inequality in (4.28) follows from Lemma 4.12 under the assignments A = 
V~X(( n ~ 2 \ C = \C\( ( - n ~ 1 \ and B\ = &(A — A)Q n- . We also use the induction assumption. 
Additionally, ( (n) (X,T) = 1 = C (n_1) (A,T); thus, C (n) < C (n_1) on G. □ 

In the next two theorems, we answer the question motivated by Theorem 4.13, that 
is, we determine the limiting value of the decreasing sequence ^P( n \ First, we show that 
lp( n ) is bounded below by Ff3, in which (3 solves 

Pt+(a- ab(X - \))(3 X + h 2 (X - A) 2 /3 AA - A/3 = 0, 
/3(A,T) = 1. 

Intuitively, the function /3 is less than (p because we replaced the square root in </?'s pde 
with the square root of the first term. Later, we show that equals Fj3 in the limit. 

In other words, lim^oo -ip^ = f3. 
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In order to prove that ip( n ) > F0, we require the following lemma. 

Lemma 4.14. The function (3 given by (4.29) is nonincreasing with respect to X. 

Proof. Differentiate /?'s equation with respect to A to get an equation for / = (3\. The 
function / solves 



ft + (a x -ab-X)f + (a- ab(X -X)+b 2 (X- X))f x + ^b 2 (X - X) 2 f xx -(3 = 0, 

/(A,T)=0. 

(4.30) 

Thus, /'s equation is of the form in (3.4) with g given by 



g(X, t, v,p) = (a x -ab- X)v + (a- ab(X - A) + b 2 (X - X))p - (3. (4.31) 

Because of the growth condition assumed for a x in Assumption 3.4, it is straight- 
forward to show that g satisfies the one-sided Lipschitz condition in (3.2) with growth 
conditions in (3.3). Define a differential operator C on Q by (3.4) with g given by (4.31). 
Because / solves (4.30), we have Cf = 0. Denote by the function that is identically 
on G; then, CO = —f3 < = Cf. In addition, /(A,T) = 0; thus, Theorem 3.2 implies that 

f = Px<0. □ 
Theorem 4.15. 

lim -P^(r, A, t) > F(r, t)(3(X, t). (4.32) 

Proof. It is enough to show that the solution f3 of (4.29) is a lower bound of -^(p( n \ We 
proceed by induction. First, show that (3 < (p = (p^ 1 '. Define a differential operator C on 
Q by (3.4) with g given by (3.8). Because <p solves (2.18), dp = 0. Also, 

£/? = a| ^Jb 2 {X - X) 2 (3\ + A/3 2 - b(X - A) \p x 1 1 > = C<p. (4.33) 

In addition, /3(A,T) = ip(X, T) = 1; thus, Theorem 3.2 and Lemma 3.3 imply that (3 <<p>. 

Next, assume that (3 < ■^i^ n ~ 1 \ and show that (3 < ■ Recall that the function 

£(»*) = solves (4.24). Define a differential operator V on Q by (3.4) with g = g n 

given in (4.27). Recall that (n - l)C (n_1) = v? (n-1) in (4.27). Then, £>C (n) = 0, and 
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Vp = p t + a(3 x + l -b 2 {\ - X) 2 (3 XX - X (n@ - ^ n ~^) 
+ a^b 2 (X-X) 2 Pl + ^X(n/3 -ip^- 1 )) 2 

= A (V n_1) - (n - + a | ^6 2 (A - A) 2 /3 2 + l\ (n/3 - v^" 1 )) 2 - b(X - A) \(3 X | J 
> = V( {n) . 

(4.34) 

Also, C (n) (^T) = (3{\T) = 1; thus, Theorem 3.2 and Lemma 4.1 imply that (3 < C (n) = 



n 



The next theorem tightens the result of Theorem 4.15 and shows that we have equality 
in (4.32). Consider the solution to the following pde: 



| 7 <»> + (a - ab(X - A)) 7 W + V(A - A) 2 7 ^ - (nX - av^x) ( 7 <»> - ^) = 0, 
\ 7 (n) (A,T)=n, 

(4.35) 

in which 7 (°) = 0. We proceed by showing in a series of lemmas that 7 < - n - ) > ip( n > for 
?i > 0. Intuitively, 7 ^ n -* is greater than because we replaced the square root in (p( n > : s 
pde with the sum of the square roots of the two terms. Finally, we show that ^ r y^ n ' — (3 
goes to zero as n goes to infinity. 

Lemma 4.16. The function 7 < - n - ) given by (4.35) is nonincreasing with respect to X and 
7 (n+i) > 7 (n) j orn > _ 

Proof. The proof of that 7 ^ < is similar to the proof that cp x n ' < in Theorems 3.7 
and 4.4, so we omit the details. Similarly, 7 ( n+1 ) > 7 ( n ) follows as in the proof of Lemma 
4.3. □ 



Lemma 4.17. ^ > v? (n) for n>0. 
Proof The result is true for n = because 

7 (o) = ^(0) 

= 0. Suppose for n > 1, we have 
^(n-i) > an( ^ s } low that 7 ( n ) > </?( n ). Define a differential operator C on Q by (3.4) 

with g = g n given in (4.3). Then, £</?( n ) = 0, and 
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£ 7 (») = ab(X - A) 7 J n) - nX ( 7 (n) - + (nX - aV^X) ( 7 (n) - 7 (n_1) ) 



+ ay b 2 (X - A) 2 ( 7 < n) ) + nA ( 7 (») - y^- 1 )) 
< (nA - av^A) (V n_1) - 7 (n_1) ) < = 



(4.36) 



In the first inequality in (4.36), we use the facts that 7 j^ < and 7 ( n ) > 7 (n_1 ), and 
we use the subadditivity of the square root function. In the second inequality, we use the 
induction assumption. We also have the terminal conditions 7^) (A, T) = n = (p( n \X,T); 
thus, Theorem 3.2 and Lemma 4.1 imply that 7^) > ip^ n \ □ 



Now, we are ready to prove the main theorem of this section. 
Theorem 4.18. 



lim -P( n) (r,A,t) = F(r,t)p(X,t). 



(4.37) 



Proof. By Theorem 4.15 and Lemma 4.17, the theorem is proved if we show that ^7^ — (3 
goes to zero as n goes to infinity because ^7^ — (3 > ^<p^ — (3 > 0. 

Define $( n ) on G by $^ n - ) = ^7^ — /3, so the theorem is proved if we show that 
lim^oo $ (n \X,t) = 0. The function ^ solves the recursion 



» 



+ (a - «6(A - A))$i n) + |& 2 (A - A) 2 < n A } - (nA - av^) $ (n) 
k $ (n) (A,T) =0, 



(4.38) 



with < <I>( 1 ' ) = 7^ 1 ) — j3 < 1. From (4.38) and the Feynman-Kac Theorem, we deduce the 
following expression for <I>( n ) in terms of $( n_1 ): 



$ (n) (A,t) = «E 

+ (n- 1)E 




— (3(X s ,s)e~Jt v»-«-«v»-«;-™ ds 
t 1 n 



A t = A 



ds 



X t = X 



(4.39) 
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in which A follows the process dX s = (a—ab(X s —X))ds+b(X s —X)dWg, with W S A = W^+as. 



The process W x is a standard Brownian motion with respect to the probability space 
(fi,jF, P), in which |p = e ~ aW T-\oP-T _ denotes expectation with respect to P. Note 
that a is analogous to the bond market's price of risk q in (2.4). 

Suppose $ (n_1) (A,t) < K n _i on G for some n > 2. Then, we can bound $ < - n - ) as 
follows: 




du 



ds 



n 



\t = \ 



+ (n-l)K n _ 1 E 



X t = X 



in which we use the fact that < (3 < 1. Define / (n ) and /i (n ) for n > 2 by 



/<">(A,f) = oE 



^ njX s e~ K^-^^Z^ds 



\t = A 



and 



(4.40) 



(4.41) 



/i (n )(A,t) = E 



A, = A 



(4.42) 



Thus, inequality (4.40) is equivalent to 



^(V) ^ -^/^(^^ + ^^n-l/i (n) (A,t). 



(4.43) 



After proving the following two lemmas that give us bounds on f^ and in (4.41) 
and (4.42), respectively, we finish the proof of Theorem 4.18. 

Lemma 4.19. The function f^ defined by (4.41) is bounded above by J = for all 

n>2. 

Proof. By the Feynman-Kac Theorem, f^ n > solves the linear pde 



fl n) + (a - ab(X - A))/<»> + \b 2 (X - A) 2 /ff - (nA - aVnJ)f^ = 
/ (n) (A,T) =0. 



-an 



Vx, 



{AAA) 
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Define a differential operator C on Q via /( n )'s equation. Let J denote the function that is 
identically equal to J. Thus, Cf^ = 0, and CJ = - (nX - aVnX^j J + anVX < = Cf^ 
for n > 2. We also have f( n \X,T) = < J; thus, from Theorem 3.2, we conclude that 
/( n ) < J on G for n > 2. □ 

Lemma 4.20. TTie function defined by (4.42) bounded above by 1 for all n > 2. 
Proof By the Feynman-Kac Theorem, ft( n ) solves the linear pde 



hi n) + (a- ab(X - A))/ii n) + \b 2 (X - A) 2 /ig - (nX - ay/nX)h^ = -(nX - av^A), 

/i (n) (A,T) = 0. 

(4.45) 

Define a differential operator £ on ^ via /i^'s equation; thus, ChS n ^ = = £1. We also 
have /i^ n - ) (A,T) = < 1; thus, from Theorem 3.2, we conclude that < 1 on G for 
n > 2. □ 

End o/ Proo/ of Theorem 4. 18. We have shown that if < K n _ l5 then $( n ) < K n , 
in which for n > 2, 

7 77 — 1 

K - = -T72 + K n-U (4-46) 

with iiTi = 1. In hindsight, we could have used induction to show that if < K n -\, 

then < K„ via a comparison argument with a differential operator based on the pde 
in (4.38). However, we have chosen to leave the proof as is because of its constructive 
nature in demonstrating the origin of (4.46). 
Define L n = nK n ; thus, 

L n = L n _! + ^=, 77 > 2, (4.47) 

V77 



from which it follows that 



i=2 

Finally, we have 



" J f n dx 

L n = l+} _ < l + J / _ < l + 2 Jv^, ?7 > 2. 
t^Vi Ji Vx 



(4.48) 



1 2 7 

& n \X,t)<K n <- + —, n >1, (4.49) 
?7 
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The right-hand side of inequality (4.49) goes to zero as n goes to oo; thus, <3>( n ) (A, t) goes to 
zero as n goes to oo. In other words, lim n ^ 00 ^p( n \\,r,t) = F(r,t)f3(\,t), as we wished 
to show. □ 

We have the following corollaries of Theorem 4.18. The first parallels the observation 
that the right-hand side in (4.20) is simply EX if b = 0. 

Corollary 4.21. // b = 0, then lim^oo ^P (n) (A,r,t) = F(r,t)^ a0 in which <p a0 is the 
physical probability of survival given in (3.15). 

Proof. From Theorem 4.18, we know that -^<p^ n ' goes to (3 as n goes to infinity. Therefore, 
the corollary follows because (3 = ip a0 when 6 = 0, which is clear from (4.29) by setting b 
equal to 0. □ 

Corollary 4.22. If b is uniformly bounded below by k > 0, then lim n ^ 00 ip( n )(A, r, t) > 
F(r, t) (fi a0 , with equality only when t = T . 

Proof. This result follows from the fact that (3 > </? a0 , with equality only when t = T. 
Indeed, define a differential operator C on Q by (3.4) and (3.8) with a = 0. Thus, C(p a0 = 0, 
and 

C/3 = ab(X - X)Px < = &p a0 . (4.50) 

Also, </? a0 (A,T) = /3(A,T) = 1; thus, Theorem 3.2 and Lemma 3.3 imply that <p a0 < (3. 

To show that </? Q ° < (3 for t G [0, T), first consider the pde of (Px°- ^ ne function 
/ = (p%° solves 



ft + (ax ~ A)/ + (a + b 2 (X - X))f x + -b 2 (X - A) 2 /aa - V> a0 = 0, 
/(A,T) = 0. 

From the linear pde in (4.51) and the Feynman-Kac Theorem, we deduce that 



(4.51) 



V f(X,t) = -B 



J ^\X s ,s)e-$t {X -- ax{X -> u))du ds 



A, = A 



(4.52) 



in which A follows the diffusion dX s = (a + b 2 (X s - X))ds + b(X s - X)dW^ with = 
Wg— Jq b(u)du. The process W x is a standard Brownian motion with respect to a suitably- 
defined probability space, and E denotes expectation on that space. Note that v 3 " ^ t) < 
for A > A and for t E [0, T) because </? a0 (A, t) > on that domain by the representation in 
(3.15). 
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Next, consider the pde of B = (3 — (p : 



1 



B t + (a - ab(X - X))B X + -b\X - X) 2 B XX - XB = ab(X - A)^f , 
B(X,T) = 0. 

From the linear pde in (4.53) and the Feynman-Kac Theorem, we deduce that 



(4.53) 



B(X,t) = -qE 



b(s)(X s -X)cpf(X s ,s)e 



■j X u du 



ds 



X, = X 



(4.54) 



in which A follows the diffusion from the paragraph following (4.39) and E is as in (4.39). 
Note that B(X,t) > for A > A and for t e [0, T) because </?^ (A, t) < on that domain 
by the representation in (4.52). □ 

After all this work, we can finally decompose the risk charge P — P a0 - first mentioned 
in the Introduction, then briefly in Section 2.2, and again following Corollary 3.9 - into 
its component risk charges: one for a finite portfolio and another for stochastic mortality. 
More generally, we decompose the per-risk risk charge \P^ — P a0 when the insurer holds 
a portfolio of n risks. Recall from Theorem 4.18 that lim n ^ 00 ^P^ = Fj3. Therefore, 
define -^P^ — F/3 as the risk charge (per risk) for holding a finite portfolio, and define 
F/3 — P a0 as the risk charge for stochastic mortality even after selling to an arbitrarily 
large group. Thus, we have 



l p (n) _ pao = ^Ip(») - F^j + (F(3 - P a0 ) = F {^P {n) -P^j+F(p- ip a0 ) , (4.55) 

in which the risk charge for stochastic mortality, namely F(f3 — <^ Q °), is zero if b = by 
Corollary 4.21 and is positive (for t < T) if b > k > by Corollary 4.22. 

5. Summary and Conclusions 

We developed a theoretical foundation for valuing mortality risk by assuming that 
the risk is "priced" via the instantaneous Sharpe ratio. Because the market for pure en- 
dowments is incomplete, one cannot assert that there is a unique price. However, we 
believe that the price that our method produces is a valid one because of the many de- 
sirable properties that it satisfies. In particular, we studied properties of the price for n 
conditionally independent and identically distributed pure endowment risks. In Theorem 
4.11, we showed that the price is subadditive with respect to n, and in Theorem 4.13, we 
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showed that the risk charge per person decreases as n increases. We also proved that if 
the hazard rate is deterministic, then the risk charge per person goes to zero as n goes to 
infinity (Theorem 4.20 and Corollary 4.21). Moreover, we proved that if the hazard rate is 
stochastic, then the risk charge person is positive as n goes to infinity, which reflects the 
fact that the mortality risk is not diversifiable in this case (Theorem 4.20 and Corollary 
4.22). Additionally, in equation (4.55), we decomposed the per-risk risk charge into the 
finite portfolio and stochastic mortality risk charges. Because of these properties, we an- 
ticipate that our pricing methodology will prove useful in pricing risks in other incomplete 
markets. 

In addition to addressing the problem of the breakdown in the law of large numbers, 
the study of such a mortality risk premium might breathe new life into the analysis of 
exotic options that are embedded within insurance and pension contracts, which have 
traditionally been viewed as being out-of-the-money and, hence, valueless. Indeed, our 
theoretical framework - which incorporates a risk measure for mortality risk - could be 
used to evaluate the risks inherent in a number of recent industry trends, namely (i) the 
emergence of longevity-linked bonds, (ii) capacity constraints in the immediate annuity 
market, and (iii) the decline of defined benefit pension plans. 

Acknowledgements 

We thank Erhan Bayraktar, Kristen S. Moore, Jeffrey Rauch, Keith Promislow, and 
Zhengfang Zhou for their valuable help. 

References 

Biffis, E. (2005), Affine processes for dynamic mortality and actuarial valuation, Insurance- 
Mathematics and Economics, to appear. 

Bjork, T. (1998), Arbitrage Theory in Continuous Time, Oxford University Press, Oxford. 

Blanchet-Scalliet, C, N. El Karoui, and L. Martellini (2005), Dynamic asset pricing theory 
with uncertain time-horizon, Journal of Economic Dynamics and Control, 29: 1737- 
1764. 

Boyle, P. P. and M. Hardy (2003), Guaranteed annuity options, AS TIN Bulletin, 33: 125- 
152. 

Brennan, M. J. and E. Schwartz (1976), The pricing of equity-linked life insurance policies 
with an asset value guarantee, Journal of Financial Economics, 3 (1): 195-213. 

Cairns, A. J. C, D. Blake and K. Dowd (2004), Pricing framework for securitization of 
mortality risk, working paper, Heriot-Watt University. 



36 



Cox, S. H. and Y. Lin (2004), Natural hedging of life and annuity mortality risks, Journal 
of Risk and Insurance, to appear. 

Dahl, M. (2004), Stochastic mortality in life insurance: Market reserves and mortality- 
linked insurance contracts, Insurance: Mathematics and Economics, 35: 113-136. 

DiLorenzo, E. and M. Sibillo (2003), Longevity risk: Measurement and application per- 
spectives, working paper, Universita degli Studi di Napoli. 

Gerber, H. U. (1979), Introduction to Mathematical Risk Theory, Huebner Foundation 
Monograph 8, Wharton School of the University of Pennsylvania, Richard D. Irwin, 
Homewood, IL. 

Gerber, H. U. and E. S. W. Shiu (1994), Option pricing by Esscher transforms (with 
discussions), Transactions of the Society of Actuaries, 46: 99-191. 

Karatzas, I. and S. E. Shreve (1991), Brownian Motion and Stochastic Calculus, second 
edition, Springer- Verlag, New York. 

Lamberton, D. and B. Lapeyre (1996), Introduction to Stochastic Calculus Applied to 
Finance, Chapman & Hall/CRC, Boca Raton, Florida. 

Lee, R. D. and L. R. Carter (1992), Modeling and forecasting U.S. mortality, Journal of 
the American Statistical Association, 87 (419): 659-671. 

Milevsky, M. A. and S. D. Promislow (2001), Mortality derivatives and the option to 
annuitize, Insurance: Mathematics and Economics, 29: 299-318. 

Norberg, R. (2004), Vasicek beyond the normal, Mathematical Finance, 14 (4): 585-604. 

Olivieri, A. (2001), Uncertainty in mortality projections: An actuarial perspective, Insur- 
ance: Mathematics and Economics, 29: 231-245. 

Protter, P. (1995), Stochastic Integration and Differential Equations, Applications in Math- 
ematics, 21, Springer- Verlag, Berlin. 

Royden, H. L. (1968), Real Analysis, second edition, Macmillan, New York. 

Schrager, D.F. (2005), Affine stochastic mortality, Insurance: Mathematics and Economics, 
to appear. 

Smith, A., I. Moran, and D. Walczak (2003), Why can financial firms charge for diversifiable 
risk?, working paper, Deloitte Touche Tohmatsu. 

Soininen, P. (1995), Stochastic variation of interest and mortality, Proceedings of the 5th 
AFIR International Colloquium, 871-904. 

Walter, W. (1970), Differential and Integral Inequalities, Springer- Verlag, New York. 



37 



Windcliff, H., J. Wang, P. A. Forsyth, and K. R. Vetzal (2005), Hedging with a correlated 
asset: Solution of a nonlinear pricing pde, working paper, University of Waterloo. 

Zariphopoulou, T., (2001), Stochastic control methods in asset pricing, Handbook of 
Stochastic Analysis and Applications, D. Kannan and V. Lakshmikantham (editors), 
Marcel Dekker, New York. 



38 



